Gap solitons and Bloch waves of interacting bosons in one-dimensional optical 
lattices: Prom the weak to the strong interaction limits 
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We study the gap solitons and nonlinear Bloch waves of interacting bosons in one-dimensional 
optical lattices, taking into account the interaction from the weak to the strong limits. It is shown 
that composition relation between the gap solitons and nonlinear Bloch waves exists for the whole 
span of the interaction strength. The linear stability analysis indicates that the gap solitons are 
stable when their energies are near the bottom of the linear Bloch band gap. By increasing the 
interaction strength, the stable gap solitons can turn into unstable. It is argued that the stable gap 
solitons can easily be formed in a weakly interacting system with energies near the bottoms of the 
lower-level linear Bloch band gaps. 
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I. INTRODUCTION 

Recent development of trapping and cooling techniques 
has enabled experimental realizations of Bose-Einstein 
condensation (BEG) of bosonic atoms and molecules in 
optical lattices When the temperature is low and 
density is high, strong interaction between particles can 
result in significant nonlinearity in these systems. Associ- 
ated with the periodicity and nonlinearity, there exist two 
important waves in these systems, namely Bloch waves 
and gap solitons. Bloch waves, which exist in both linear 
and nonlinear periodic systems, are extensive and spread 
over the whole space On the contrary, gap solitons, 
which are spatially localized atomic wave packets, exist 
only in a nonlinear periodic system Q. In particular, 
a class of solitons called the fundamental gap solitons 
(FGSs), have the major peak well localized within a unit 
cell 13]. The solitons with two peaks of opposite signs 
within a unit cell arc called the subfundamental solitons 
0. The relationship between the nonlinear Bloch waves 
(NLBWs) and the gap solitons is a topic of considerable 
interest. 

Both FGSs and NLBWs can be simultaneously ob- 
tained by solving numerically the nonlinear Schrodinger 
equation of the system. For the weakly interacting 
one-dimensional (ID) periodic Bosc system, it has been 
shown that NLBW can be regarded as the superposition 
of FGSs in an infinitely-long chain H, Q- This, so-called 
the "composition relation", leads to the prediction that 
there are n families of FGSs in the nth band gap of the 
corresponding linear periodic system. It also implies that 
a class of solutions similar to the Bloch waves can be built 
from it. It is interesting to see whether the composition 
relation remains correct in a strongly interacting system? 
To what extent the FGSs and NLBWs change when the 
interaction changes? 

As is well known, mean-field theory typically does not 
work well for a ID system, except in the very weakly 
interacting regime. The enhanced quantum fluctuation 
is significant in ID quantum systems (^-[Hl which ex- 



hibit fascinating phenomena significantly different from 
their three-dimensional counterparts. On the other hand, 
for a single-component Bose gas in a harmonic trap, 
it has been shown that with the increase of the repul- 
sive interaction, the density profile evolves continuously 
from a Gaussian-like distribution of bosons to a shell- 
structured distribution of fcrmions, called the Tonks- 
Girardeau (TG) gas [13, [13 ■ Thus when the interac- 
tion is strong, non-perturbative methods such as the 
Bose-Fermi mapping [TJ] or the Bethe ansatz [l^ [l^ 
need to be used to characterize the features of the sys- 
tem properly. Recently ID harmonic trapped spinless 
Bose systems with a repulsive J-function interaction have 
been solved in the hmit of -)• oo [13, [13, Ell • More- 
over, considering the two-component Bose gas with spin- 
independent interactions in the absence of the external 
potential, the ground-state energy density function has 
been extracted from the Bethe-ansatz solution. With the 
local-density approximation, a modified coupled nonlin- 
ear Schrodinger equations for the ground state has been 
obtained [l9j . Using this kind of nonlinear Schrodinger 
equation, one is able to study the interplay between the 
periodicity and the nonlinearity in a ID Bose system, es- 
pecially for the interaction from the weak to the strong 
limits. 

In this paper, we attempt to study the NLBWs and 
FGSs of an interacting Bose system in a ID optical lat- 
tices. In particular, we are interested in how the NLBWs 
and FGSs of the system change when the interaction is 
changed from the weak to the strong limits, and whether 
the composition relation remains correct upon the change 
of the interaction. It will be shown that the composition 
relation remains valid in all interaction regimes. How- 
ever, linear stability analysis indicates that stability of 
the soliton waves will change with the change of the non- 
linearity. 

The paper is organized as follows. In Sec. |TI1 we intro- 
duce the model equation for a ID interacting Bose system 
in optical lattices to which the interaction can be changed 
from the weak to the strong limits continuously. Sec. IIIII 
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is devoted to compute the ground-state wave functions 
and energies of the system for the interaction from the 
weak to the strong hmits. It is aimed to give a direct 
insight of the interaction dependent ground-state prop- 
erties. By giving as much as possible evidence numeri- 
cally in Sec. |IVl we show that the composition relation 
exists between the NLBWs and FGSs in all interaction 
regimes. A generalized composition relation between the 
high-order solitons and multiple periodic waves is also 
shown. In Sec. |V1 the stabilities of various soliton waves 
are investigated upon the changes of the interaction and 
the strength of the periodic potential. Sec. I VII is a brief 
summary. 



Bloch waves and gap solitons in optical lattices all have 
the form '^{x,t) = $(a;) exp(— i/ii). Substituting it into 
the time-dependent nonlinear Schrodinger equation P]). 
one obtains the following dimensionless time- independent 
nonlinear Schrodinger equation 



V cos(a;) + F (p) 



$(a;) ^ ^l'^>{x). (5) 



Eq. ([5]) is the starting point of the calculations through- 
out this paper. 



III. GROUND-STATE PROPERTIES 



II. MODEL EQUATION 

We consider a ID periodic Bose system described 
by the following time-dependent nonlinear Schrodinger 
equation 



' dt 



F{p) 



^{x,t), (1) 



where V{x) = ucos(^a;) is the periodic potential with 
A the lattice constant and v the strength and F (p) is 
responsible for the interaction energy. The particle num- 
ber is determined via N — pdx with p = \'^\^ ■ To go 
beyond the weak interaction limit, we use the following 
nonlinear form for F {p) [l9j . 



F{p) 



2^2 



2m 



[3e (7) - -fd^e (7)] , 



where 



e(7) = 



q"f + P7 



(2) 



(3) 



The results of Eqs. ^ and ([3]) were recently obtained by 
Hao and Chen based on Bethe Ansatz [l^ . They can be 
equivalently obtained by solving exactly the equation for 
6(7) [see, for example, Eq. (6) in Ref. [13], as done, for 
instance in Ref. Here ^ = c/ p with c = mg/h^ 

[g is the scattering length) and p = —5.0489470 and 
q = -20.8604983 are the fitting parameters In the 
weakly (7 ^ 1) and strongly (7 1) interacting limits, 
one has the asymptotic forms 



F{p) 



mg 
2m ^ 



7 < 1 
7 > 1. 



(4) 



For convenience, dimensionless scaling will be made for 
the length and the energy @,0|- Position x is to be scaled 
in the unit of A/(27r); Periodic potential V{x), interaction 
energy F{p), and the chemical potential p, are scaled in 
the unit of SEr with Er = fPn^ /2mA'^ being the recoil 
energy and m the atom mass. Wave functions of both 



By solving Eq. ([5]) numerically, one is able to obtain 
the ground-state density profile p{x) and the correspond- 
ing chemical potential p. It is remarked that to solve 
Eq. (IHl), we first differentiate it using the finite-element 
method along with the periodic boundary condition [20| , 
and then evaluate several hundreds of steps in imaginary 
time until the lowest chemical potential p is reached. In 
the calculation, the system is taken to be 16 lattice con- 
stant long {x ranges from — IGtt to IGtt), the strength 
of periodic potential v = 1.5, and the particle number 

— 48 (average 3 particles per unit cell) . Figure [TJa) 
shows the density profile p{x) within a unit cell. In view 
of Fig. [U^a), similar to the case of the Bosc gas in a single 
harmonic trap, p{x) evolves from the Bose distribution 
to the Fermi-like one when the interaction constant c is 
increased. In the weakly interacting regime (c = 0.2), the 
density profile displays a Gaussian-like Bose distribution. 
When c is increased (see, for example, the case of c = 2), 
density profile decreases at the center while increases at 
the two sides. In the extremely strong interaction regime 
(c = 500), due to the strong enhancement of particle 
tunnelling rate across neighboring wells, density profile 
at the boundary increases significantly. In such limit, so- 
called the TG gas, the system will behave similar to the 
noninteracting fermions f2l| . 

In order to compare the results obtained by the 
imaginary-time method, we also calculate the (exact) 
density profiles for the TG gas using the Bose-Fermi map- 
ping method. In this case, Eq. ^ is first solved by the 
finite-element method by setting F{p) = 0. The density 

profile is then calculated via ptg{x) — J2n=i 
with $„(x) denoting the n-th energy eigenstate of the 
Hamiltonian. It turns out that the curve of ptg^x) 
agrees well with the p{x) obtained by the imaginary-time 
method for the large c = 500 case [see Fig.[T]Ja)]. We no- 
tice that in addition to the central peak, pxG (2;) exhibits 
two humps at two sides near the center. Since the bot- 
tom of a single well behaves as an harmonic potential, 
the density ptg{x) obtained by us is very similar to the 
Tonks density for three bosons in an harmonic potential. 
Three is exactly the number of particles put on average 
in a single well. It has been shown that by increasing 
the number of particles, the number of humps increases 
correspondingly while their amplitude decreases (22j . 
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FIG. 1. (Color online) Panel (a): The interaction constant c- 
dependent ground-state density distribution p{x) within one 
unit cell. Dimensionless x ranges from to 2tt for one unit 
length. The large c = 500 case is seen to behave like the non- 
interacting fermions (blue line, see text). Panel (b): The 
c-dependent ground-state chemical potentials calculated in 
the weak limit {F{p) ~ blue line), in the strong limit 

{F{p) ~ 1$!'', red line), and for a full F{p) (circle line) re- 
spectively. The inset shows an enlarged view of the curves 
at small c's. The mark corresponding to c = 0.0009, is 

what used in Ref. 0| to discuss the nonlinear bands. 



It is well known that the dilute Bose-condensed gas 
can be described by the Gross-Pitaevskii (GP) mean-field 
(MF) theory [2^. In fact, MF approximation is only ap- 
propriate for a long-wavelength theory. For short-range 
repulsive interactions, MF theory fails in dimension d < 2 
to which the nonlinear term F{p) ~ in the GP MF 
theory should be replaced by F(p) ^ |<I>|^ from the renor- 
malization group analysis |24| . One may ask to what ex- 
tent the current system under studies will approach the 
above two (weak and strong) limits upon the change of 
the interaction constant c. To uncover it, in Fig. [Ijb) 
we study the ground-state chemical potential fi as the 
function of interaction constant c. In view of Fig. [ijb), 
one sees that the curve of F{p) ^ |$p matches well the 
one using full F{p) when c < 1. Therefore the present 
Bose gas can be taken to be in the weakly interacting 
limit when c < 1. It is noted that the point c = 0.0009, 
marked by "+" sign, is what used in Ref. to discuss the 
nonlinear bands. The difference between the two curves 
increases with the increase of c. On the other side, the 
curve using full F{p) matches well the one of F{p) ^ |<I>|^ 
when c > 4. Therefore the present Bose gas can be con- 
sidered to be in the strong TG gas limit when c > 4. 



IV. GAP SOLITONS AND BLOCH WAVES 

When the nonlinear interaction term F(p) is set to 
zero, Eq. J5]) turns into the well-known linear Mathieu 
equation |^. Its eigenfunctions are linear Bloch waves 
(LBWs) and eigenvalues form the linear Bloch band 
(LBB). Physical solutions are forbidden in the band gaps 
between neighboring bands. However, when the nonlin- 
ear interaction exists [F{p) ^ 0], two kinds of solutions 
can exist in the LBB gaps. One is the gap soliton and 



another is the NLBW. Both of them can share the same 
chemical potential in the LBB gaps. It has been illus- 
trated in Ref. Q that for the weak hmit, F{p) ^ 
the first nonlinear Bloch band (NLBB) can be lifted into 
the second, third, and higher LBB gaps with the increase 
of the scattering length g. Similarly the second NLBB 
can also be lifted into the third and higher LBB gaps. 
Consequently in the weak limit, NLBB in the nth LBB 
gap should include n-1 branches which are lifted from the 
lower levels. There is only one branch of NLBB that can 
develop in the lowest LBB gap 0. The present issue is 
to check whether the above results remain correct in a 
more realistic model such that the nonlinear interaction 
term is described by the full F{p). Our following numer- 
ical results will provide strong evidences that the above 
results remain valid for the whole span of the interaction 
regimes. 

It is worth noting that as shown in Ref. [2^ , mean-field 
theory of Eq. ([T]) may not correctly describe the physics 
of a ID impenetrable Tonks gas. This means that in 
the Tonks limit, a solitonic solution may not really exist. 
While this poses a limit on the solitonic solutions that 
we have obtained, it is still interesting to experimentally 
verify our theoretical predictions. 



A. Interaction-dependent gap solitons and 
nonlinear Bloch bands 

We first set -F(p) = to solve the linear Schrodinger 
equation ([5]) exactly by the finite-element method. 
The periodic potential strength v-dependcnt lowest four 
bands and band gaps are shown in Fig. [2ja). As shown, 
with the increase of f, the continuum energy spectrum 
turns into energy bands with gaps. When v is large 
enough, some (lower) bands reduce to highly-degenerate 
thin levels. We next retain F{p) to solve the nonlin- 
ear Schrodinger equation ^ numerically. In particular, 
we arc interested in the solutions of FGSs (or NLBWs) 
with the energy p falling into the LBB gap regions. 
The numerical solutions arc obtained by differentiating 
Eq. ([5]) on a finite difference grid to obtain a coupled alge- 
braic equations and solve it with the Newton-relaxation 
method 

Figure[2l^b) shows particle number of gap solitons as 
the function of p for four different interaction constant 
c's. As shown, for example, when c = 3 close to the 
strong limit (solid lines), the first NLBB develops from 
the first LBB and is lifted with the increase of N . The 
so-called NLBB lifting is simply due to the fact that the 
larger N is, the larger the nonlinearity and hence the 
corresponding p are. The first NLBB enters the second 
LBB gap at iV ~ 2.38 and when N ~ 4.92 it further 
enters the third LBB gap. For the second NLBB, it de- 
velops from the second LBB at A'' ~ 0.14 and is lifted 
in a similar manner with the increase oi N. It enters 
the third LBB gap when it crosses the critical value of 
N ~ 3.06. It turns out that the third NLBB only dcvel- 
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FIG. 2. (Color online) Panel (a): Energies /i of the lowest four 
linear Bloch bands and band gaps versus the strength of the 
periodic potential v. fi and v are in units of 8Er with Er the 
recoil energy (see text). The dashed line denotes u = 2 used 
in the calculations of panel (b). Panel (b): Particle number A*' 
of FGSs as the function of chemical potential jj, for different 
interaction constant c — 1.7, 3, 50, and 500, respectively. The 
points marked by red "+"sign are those studied in Fig. [4] 
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FIG. 3. (Color online) Amplitudes of the first-family FGS as 
the function of the interaction constant c. For all curves, par- 
ticle number is fixed at A'' = 3 and periodic potential strength 
is n = 2. For comparison, the case of the strong interaction 
limit, F{p) ~ I*!"!*, is also shown. 



smaller. 



ops a.t N > 0.66, slightly above the third LBB. 

For relatively weaker interaction case, c = 1.7 (dotted 
lines), the NLBB is seen to correspond to larger N for the 
same fi, as compared to those of the c = 3 case. This is 
because for the smaller interaction constant c, one needs 
a larger N to achieve the same nonlinear effect. We have 
also calculated the extremely strong interaction cases (for 
c = 50 and 500) for energies within the second LBB gap. 
The lines of the two cases are seen to coincide and indi- 
cate that the current Bose system is reduced to the TG 
gas when the interaction constant c is large enough. 



B. Interaction-dependent amplitude of gap solitons 

The soliton wave occurs only in the nonlinear system 
and has no classical correspondence. One may guess that 
the amplitude of soliton wave increases with the increase 
of the nonlinearity. However, the truth is in the oppo- 
site direction. As shown in Fig. [3l numerical results of 
the first-family FGS waves show that their amplitudes 
actually decrease as c increases. For all curves in Fig. [3l 
particle number is fixed at iV = 3 and periodic potential 
strength is v = 2. For comparison, we have also shown 
the case for the strong interaction limit [F{p) ~ l^"!^]- 

In fact, FGSs originate from the nonlinearity and peri- 
odicity of the system. On one hand, FGSs can be taken 
roughly as the discrete eigenstatcs of the system in a sin- 
gle potential well. On the other hand, FGSs are eigen- 
statcs of the Schrodinger equation ([5]) if [V{x) + F{p)] 
is taken approximately as the effective field of a single 
well. When the interaction constant c is increased and 
for the particle number N fixed, the confinement due to 
the effective single well will become weaker. This means 
that the wave function will become more dispersed and 
correspondingly the amplitude of the FGSs will become 



C. Composition relation at strong interactions 

As shown previously, NLBB can be viewed as the lifted 
LBB by increasing the nonlinear interaction. While LBB 
can be viewed as the evolution from the discrete energy 
levels of an individual well. Therefore it is not difficult 
to sec that NLBW belonging to the nth NLBB should 
have n-1 nodes (in the sense of the nth bound state) 
in an individual well of the periodic potential. It has 
been pointed out and proven numerically for the weak 
interaction limit [F{p) ~ that FGSs which can only 
exist in a nonlinear system are the building blocks of 
the NLBWs 0. Thus it is expected that FGS should 
behave like the bound states in an individual well of the 
periodic potential. The current issue is to verify whether 
the above conclusion remains correct when the nonlinear 
term is replaced by the full F{p). 

In Fig. m we solve and plot both NLBW (blue dashed 
line) and FGS (red solid line) for the relatively strong 
interaction (c = 3) case. It should be emphasized that 
the conclusion drawn here remains valid for much more 
stronger interaction cases. We consider both NLBWs and 
FGSs in the first, second, and third LBB gap respectively 
and for the representative we take the k points at both 
the BZ center {k = 0) and boundary (fc = 0.5). In view 
of Fig. |4i;b)&(c) with ^ = -0.65 and N = 1.30, an al- 
most perfect match is found between the NLBW and the 
FGS within one unit cell in the first LBB gap. The good 
match occurs regardless of at the center or at the edge 
of the BZ. It thus gives a strong evidence that FGSs can 
be considered as the building blocks of the NLBWs. It is 
important to note that the waves shown in Fig. |3{b)&(c) 
are belonging to the first-family FGS and hence of no 
node. In Fig. me)&(f), a good match is also found be- 
tween the NLBW and the FGS within one unit cell in 
the second LBB gap. These waves with = 0.5 and 
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FIG. 4. (Color online) Left column: panel (a),(d), and (h) 
show the first, second, and third LBBs and band gaps. Thick 
black dots denote the (fe, fi) points studied in the right column. 
Right column: panel (b) and (c) are the FGS (red solid line) 
and NLBW (blue dashed line) in the first LBB gap with A'' = 
i.30 and (fc,/.i) = (0,-0.65) (at BZ center) and (0.5,-0.65) 
(at BZ edge) respectively. Panel (e) and (f) are the FGS and 
NLBW in the second LBB gap with iV = i.52 and (fc,/i) = 
(0, 0.5) and (0.5, 0.5) respectively; panel (i) and (j) are the 
FGS and NLBW in the third LBB gap with = f .79 and 
(fc, /i) = (0, i.5) and (0.5,1.5) respectively. For all panels, 
the interaction constant is c = 3 and the periodic potential 
strength is w = 2. The (A'^, /i) points under study are also 
marked in Fig. [H^b). 
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FIG. 5. (Color online) Panel (a)-(c): Illustration of the high- 
order gap solitons (red solid line) and NLBWs (blue dashed 
line) for c = 3 and jj, — 1.63. (a) The upper power branch 
of five-peak gap solitons; (b) The lower power branch of five- 
peak gap solitons; (c) The lower power branch of eight-peak 
gap solitons. Panel (d)-(g): Illustration of the high-order 
gap solitons (red solid line) and multiple periodic waves (blue 
dashed line) for c = 0.3 and fi = —1.07. (d) Two-peak solitons 
with A'' = 6; (e) Solitons of two peaks up and two peaks down 
with A'' = 12; (f) Solutions of three peaks up and three peaks 
down with A^ — 18; (g) Solitons of four peaks up and four 
peaks down with A'^ — 24. The periodic potential strength is 
V — 2 for all panels. 



N = 1.52 arc belonging to the second-family FGS and 
hence of one node. There is a key difference between the 
NLBWs in Fig. Hfb) and (c) such that the waves are os- 
cillating in-phase (k = 0) between the neighboring wells 
in case (b), while it is out-of-phase {k = 0.5) in case (c). 
Similar consequence is also seen between Fig. a-nd 
(f) and also between Fig. IJlJi) and (j). 

As a test, we actually try to build the possible NLBWs 
using the first obtained third-family FGS in the third 
LBB gap [see Fig.|4i;i) and (j)]. These third-family FGSs 
with /i = 1.5 and N = 1.79 are obtained by the Newton- 
relaxation method with proper initial condition. To cap- 
ture the out-of-phase motion of the NLBWs at the BZ 
edge [case (j)], we have added an alternative negative sign 
to the building block. In view of Fig. HJi) and (j), the 
match between the NLBW and FGS is good except FGS 
exhibits a short tail across the cell boundary. Neverthe- 
less, one can still conclude that the composition relation 
between NLBWs and FGSs is valid to a good approxi- 
mation. 



D. Generalized composition relation 

In addition to the FGSs and NLBWs studied previ- 
ously, there are two types of waves which are also com- 
mon in a nonlinear system. One is called the high-order 
gap solitons which are gap soliton waves of multiple peaks 
over multiple unit cells |2a]. As we will show later, high- 
order gap solitons can be viewed as the truncated NLBW 
and hence NLBWs can be built from them. Another is 
called the multiple periodic waves which are defined as 
$(x) = exp(ikx)'ipkix) with ij^kix) = ipkix + 2mr) and 
n being a positive integer [27[. In the weakly interact- 
ing limit, it has been shown that composition relation 
between the FGSs and the NLBWs can be generalized 
to construct multiple periodic waves from the high-order 
gap solitons 01 • The present issue is again to see whether 
the generalized composition relation remains valid in the 
strong interaction limit. 

We have first solved the time-dependent nonlin- 
ear Schrodinger equation ([T]) numerically using the 
imaginary-time method to obtain NLBW function $(a;) 
and chemical potential fi with a given interaction con- 
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stant c. The Newton-relaxation method is then used to 
solve the time-independent nonlinear Schrodinger equa- 
tion ([5]) with the pre-obtained and the proper ini- 
tial condition. Fig. [Ha)-(c) show the five- and eight- 
peak high-order gap solitons with the interaction con- 
stant c — S. An almost perfect match is found between 
the high-order gap solitons and the NLBWs. Thus high- 
order gap solitons can also be viewed as the truncated 
NLBWs for the strong interactions. 

Fig. [5jd)-(g) show the composition relation between 
the high-order gap solitons and multiple periodic waves. 
Here c = 0.3 is used. Shown in Fig. [S{d) is the building 
block of two-peak high-order gap solitons which gener- 
ates the two-peak up and two-peak down multiple peri- 
odic wave. Fig. El^e) shows the same two-peak up and 
two-peak down multiple periodic wave that can also be 
built by the two-peak up and two-peak down high-order 
gap solitons. We further show the three-peak up and 
three-peak down and four-peak up and four-peak down 
multiple periodic waves in Fig. [S{f ) and (g) built by the 
three-peak up and three-peak down and four-peak up 
and four-peak down high-order gap solitons respectively. 
The matches between the higher-order gap solitons and 
multiple periodic waves are all seen to be good. This 
stands that for the whole range of the interactions, FGSs 
arc the most basic building blocks which can be used to 
build high-order solitons and multiple periodic waves. 



V. STABILITY OF THE SOLITON WAVES 

Stability of the FGS is another important issue when 
the interaction of the system is changed from the weak 
to the strong limit. We shall study the linear stability 
of FGS solution following the standard procedure. Since 
the unstable solution is sensitive to a small perturbation, 
one can add a small perturbation A^(x,t) to a known 
solution 

^'(a;, t) = [^{x) + A$(a;, t)] exp(-z/xt), 

where A<f>(x,t) = u(x) c:xjp{iat) + w* (x) cyij>{~ia*t). In- 
serting the perturbation into Eq. ([1]) and dropping the 
higher-order terms in one then obtains the linear 

eigenequation 



0.15 




(6) 



where C=^£ 



V{x)-F{p)-^\^{xr+f^.Une^i- 
stability of a soliton is determined by the energy spec- 
trum of the linear eigenequation Among all eigen- 
values a obtained, if there exists one of a finite imaginary 
part, the solution of $(a;) would be unstable. Otherwise, 
the solution of 4>(x) is stable. 

The stability of FGSs is investigated in Fig. [SI For the 
weak interaction case (c = 1), the linear stability anal- 
ysis indicates that the first family FGSs which develop 
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FIG. 6. (Color online) Studies of the stability of the FGSs 
in the first, second, and third LBB gap for weak (c = f) and 
strong (c = 50) interaction limits. The periodic potential 
strength used is w = 2. 



from the first LBB are stable when their chemical poten- 
tials n are near the bottom of the first LBB gap. They 
will become unstable when ^ becomes higher within the 
first band gap (with the increase of N), and enters into 
the second and third band gaps (not shown in Fig. [6]). 
Similarly for the second family FGSs which develop from 
the second LBB, they are also stable when the chemical 
potentials are near the bottom of the second band gap. 
While for the third family FGSs, they are stable only 
when fj, is extremely close to the bottom of the third 
band gap. In contract, for the strong interaction case 
(c = 50), studies of linear stability indicate that stable 
FGSs of the first and second families exist only in a much 
narrower regime with fi near the bottom of the band gaps; 
while it seems that the third-family FGSs can no longer 
be stable. Bearing the above results, it is argued that 
stable FGSs can be easily formed in a weakly interacting 
system with chemical potential near the bottom of the 
lower-level band gaps. 

For completeness, we have also investigated the sta- 
bility of two kinds of five-peak gap solitons. One is the 
upper power branch shown in Fig. Elja) and another is 
the lower power one shown in Fig. [Sljb) [l!]. The linear 
stability analysis was carried out in the first LBB gap 
for three chemical potentials, ^ = —1.3, —1.22, and —1.2 
respectively and the results are shown in Fig. [71 For 
the case oi fi = —1.3, five-peak gap solitons are stable 
for both upper [see Fig. WIS)] and lower power branches 
[see Fig. [Zje)]. This is consistent with the case of FGSs 
which are stable with fi near the bottom of the LBB 
gap. For the case of /i = —1.22 near the middle of the 
first band gap, the lower power branch of five-peak gap 
solitons are still stable [see Fig. [Tljc)], while the upper 
power branch of five-peak gap solitons become unstable 
[see Fig. [7]Jd)]. This is consistent with the one analysis 
studied in Ref. [2^. When fj, is increased to be —1.2, 
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FIG. 7. (Color online) Linear-stability studies of the five- 
peak gap solitons in the first LBB gap. Panel (a), (c), and 
(e) correspond to the lower power ones shown in Fig. [5][b), 
while panel (b), (d), and (f) correspond to the upper power 
ones shown in Fig. [Sla). Chemical potential /i = — f .3 for 
(e) and (f), -1.22 for (c) and (d), and -1.2 for (a) and (b). 
Interaction constant c = 1 and periodic potential strength is 
v = 2. 



five-peak gap solitons are unstable for both the lower 
[see Fig. [Tlja)] and upper power branches [see Fig. [Tijb)]. 
Therefore, the stable multi-peak solitons can also be eas- 
ily formed with n near the bottom of the lower band 
gaps. 

The stability of the multi-peak solitons upon the 
change of the periodic potential strength (v) is also stud- 
ied. It is found that all the unstable multi-peak solitons 
studied in Fig. [7| will eventually become stable when v 
is increased from the present value 2 to 10. The reasons 
behind it are easily understood. When v is increased, 
interaction is reduced relatively. Consequently chemical 
potential will approach relatively closer to the bottom of 



the LLB gap and gap solitons will become relatively more 
stable. 



VI. SUMMARY 

In summary, we have investigated the composition re- 
lation between the fundamental gap solitons (FGSs) and 
nonlinear Bloch waves (NLBWs) of interacting bosons in 
one-dimensional optical lattices. The main focus is to 
consider the interaction from the weak to the strong lim- 
its. Our numerical results verify that the composition re- 
lation remains correct in the whole span of the interaction 
strength. FGSs are thus the fundamental building blocks 
to build all other stationary solutions, including NLBWs, 
high-order solitons, and multiple periodic waves, in one- 
dimensional nonlinear periodic systems. By the linear 
stability analysis, it is found that the stable FGSs exist 
near the bottom of the linear band gap. By increasing 
the interaction strength, the stable gap solitons will be- 
come unstable. However, one can restore the stability of 
the unstable FGSs by increasing the strength of the pe- 
riodic potential. It is argued that the stable gap solitons 
can be easily formed in a weakly interacting system with 
the chemical potential near the bottom of the lower-level 
band gaps. 
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